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Abstract: We study the correspondence between scattering amplitudes and Wilson 
loops in three-dimensional Chern-Simons matter theories. In particular, using M = 
2 superspace formalism, we compute at one loop the whole spectrum of four-point 
superamplitudes for generic Af > 2 supersymmetric theories and find a vanishing result 
for M = 6 AB J(M) and M = 8 BLG models. This restricts the possible range of theories 
for which Wilson loops/scattering amplitudes duality might work. At two loops, we 
present the computation of the four-point ABJ scattering amplitude for external chiral 
superfields. Extending the known result for the ABJM Wilson loop to the ABJ case we 
find perfect agreement. We also discuss the dual conformal invariance of our results and 
the relationship between the Feynman diagram computation and unitarity methods. 
While for the ABJM theory dual conformally invariant integrals exactly reproduce the 
amplitude, for the ABJ case this happens only up to a residual constant depending 
on the parity-violating parameter. Finally we propose a BDS-like exponentiation for 
the amplitude based on an analogy with the four dimensional M = 4 SYM case, and 
discuss its strong coupling dual counterpart. 
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1. Introduction 

Recently, new interest has been devoted to the study of the S-matrix for the non- 
trivial sector of three dimensional Chern-Simons-matter theories which allow for a 
string theory dual description. These are the well-known J\f = 6 superconformal AB JM 
model [1] for U(N) K x U(N)_ K gauge group and the more general ABJ model [2] for 
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U(M)k x U{N)^Ki where K is the Chern-Simons level. In the large M,N limit their 
strong dual description is given in terms of M-theory on AdS 4 x S 7 /Z K background 
and, for K < N < K 5 , by a type IIA string theory on AdS 4 x CP 3 . 

The main motivation is to understand whether these theories, even if distinguished 
in nature being them non-maximally supersymmetric, share fundamental properties of 
the four dimensional M = 4 SYM theory, like integrability [3] , Yangian symmetry [4] of 
the planar physical sector and (scattering amplitudes) / (Wilson loop (WL))/ (correlation 
functions) dualities [8]— [13], [14]-[17]. While going deep into the nature of three di- 
mensional theories, the investigation of these properties should help to understand their 
actual origin and the role of the AdS/CFT in their determination. 

For the ABJM model, preliminary results can be already found in literature, con- 
cerning integrability [18]-[33], the related Yangian symmetry [34, 35] and dualities. 

At classical level, scattering amplitudes have been shown to be invariant under dual 
superconformal symmetry [36, 37] whose generators are the level-one generators of a 
Yangian symmetry [34]. At strong coupling this symmetry should rely on self-duality 
properties of type IIA string on AdS 4 x CP 3 under a suitable combination of bosonic 
and fermionic T-dualities [38, 39], even if the situation is complicated by the emergence 
of singularities in the fermionic T-transformations [40]- [44]. Considerable progress in 
this direction has been recently done in [45]. 

At quantum level, first evidence of the existence of dualities and the persistence of 
dual conformal invariance comes from recent findings on scattering amplitudes, light- 
like WL and correlators of BPS operators at one and two loops. 

At one loop both the four-point amplitude [46] and the light-like four-polygon WL 
[47, 48] vanish. For J\f = 2, 3, 6, 8 Chern-Simons-matter theories, correlation functions 
of 2n BPS operators have been computed [48] 1 . It has been proved that the one-loop 
result divided by the corresponding tree level expression coincides with the one-loop 
light-like 2n-polygon WL [47]. The identification is at the level of the integrands, 
independently of the fact that both of them eventually vanish. 

Less trivial evidence for a (scattering amplitude)/WL duality arises at two loops 
where these quantities are not supposed to vanish. Very recently, it has been proved 
that for the ABJM theory, at this order the four-point scattering amplitude [51, 52] 
divided by its tree-level counterpart coincides with the second order expansion of a 
light-like four-polygon Wilson loop [47]. 

lr The result of [48] holds also for the TV = 8 BLG theory [49, 50] described by SU(2) K x SU(2)- K 
Chern-Simons-matter theory in the large K limit. Enhancement to maximal supersymmetry could 
also be obtained for Chern-Simons levels K = 1,2. However, these values are out of the perturbative 
regime and will not be considered. 
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In this paper we give details of the calculation for the four-point scattering ampli- 
tude up to two loops and extend the result to the more general ABJ theory. 

Using Af = 2 superspace description and a direct Feynman diagram approach, 
at one loop we compute the whole spectrum of four-point superamplitudes for Af = 
2, 3, 6, 8 Chern-Simons matter conformal field theories with U(M) x U(N) gauge group. 
The result is generically different from zero, except for the Af = 6 ABJ(M) and Af = 8 
BLG cases where they all vanish. Suitably generalizing the definition of Wilson loop 
to the ABJ theory, we easily argue that it also vanishes at one loop. Therefore, we 
conclude that a (scattering amplitude) /WL duality may work only for the Af > 6 case, 
independently of the fact that conformal symmetry is present in all the theories we 
analyze. 

Focusing on ABJ models, still using a direct Feynman diagram approach, we eval- 
uate the two-loop planar scattering superamplitude of four chiral superfields, two of 
them in the bifundamental and two in the antibifundamental representation of the 
U{M) x U(N) gauge group. The result for the ratio Mf = Af loops) /At ee is 



Mf = XX 



(2e) 2 (2e) 2 2 



In 2 (|) + CU(M, N) + (9(e) 



where A = M/K, X = N/K, \jf is the mass scale and C^(M, N) is a constant depending 
on the ranks of the groups. For M = N we are back to the result for the ABJM theory 
[51, 52]. 

This result has a number of remarkable properties. 

First of all, up to an additive, scheme-dependent constant, this expression matches 
exactly the second order expansion of a ABJ light-like four-polygon Wilson loop, once 
the IR regularization is formally identified with the UV one and the particle momenta 
are expressed in terms of dual coordinates (s = x\ z and t = x| 4 ). Therefore, at least 
for the four-point amplitude, there is evidence that the following identity 

\nM 4 = ln(W 4 ) + const. (1.2) 

should hold order by order in the perturbative expansion of the two objects. 

The result (1.1) can be identified with the first order expansion of a BDS-like 
ansatz for the ABJ(M) model 



M 4 



; W«+^Ca^(ln»(f) + ^) + C(A,A) (L3) 



where C(A,A) is a scheme-dependent constant. This ansatz is exactly the BDS ansatz 
[53, 54] for Af = 4 SYM where the four dimensional scaling function has been substi- 
tuted by the three dimensional one, fcs(X, X), which is an obvious generalization of 
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the ABJM scaling function obtained from the conjectured asymptotic Bethe equations 
[19]. 

In the Af = 4 SYM case, the BDS exponentiation of scattering amplitudes holds 
also at strong coupling where, according to the Alday-Maldacena prescription [38], the 
amplitudes are given by the exponential of a minimal-area surface in the AdS 5 dual 
background ending on a light-like polygon, whose edges are determined by the particle 
momenta. In particular, since this prescription is equivalent to computing a light-like 
WL in AdS, it supports the amplitude/WL duality at strong coupling, in agreement 
with the findings at weak coupling. 

The natural question which arises is whether a similar prescription at strong cou- 
pling can be formulated for the ABJ(M) models. Motivated by the evidence in favor 
of the amplitudes/WL duality at weak coupling and BDS exponentiation, we expect 
it to be the case. In fact, focusing on the ABJM theory in the intermediate regime 
K ^ N <^ K 5 , we discuss the generalization of the Alday-Maldacena prescription to 
AdS4 x CP3. We find that, apart from a rigorous prescription for the regularization 
procedure in AdS that we have not formulated properly, the five dimensional solution 
can be adapted to the four dimensional case and the output is an expression for the 
four-point amplitude given by eq. (1.3) where the scaling function assumes its leading 
value at strong coupling, /cs(A) ~ V2X- 

For M = N, our result (1.1) coincides with the one in [51] obtained by making 
the ansatz that dual conformal invariance should hold also at loop level. Therefore our 
calculation supports that ansatz and provides a direct proof of the assumption that 
dual conformal invariance should be the correct symmetry principle to select the scalar 
master integrals contributing to the on-shell sector of the theory. 

In fact, for the ABJM case, following [51] we can rewrite the result (1.1) as a linear 
combination of scalar momentum integrals which are dual to three dimensional true 
conformally invariant integrals, well defined off-shell. As a consequence, the four-point 
amplitude satisfies anomalous Ward identities associated to dual conformal transforma- 
tions [8], as dual conformal invariance is broken in the on-shell limit by the appearance 
of IR divergences which require introducing a mass regulator. 

For the ABJ model the situation is slightly complicated by the appearance of 
a non-trivial dependence on the parity-violating parameter a = (M — N) / V MN 
in the mass-scale and in the constant C4. In fact, for M ^ N the two-loop ratio 
(1.1) can be still written as a linear combination of dual invariant integrals only up 
to an additive constant proportional to a 2 . Therefore, for the ABJ theory the dual 
conformal invariance principle combined with the unitarity cuts method is not sufficient 
to uniquely fix the amplitude already in the case of four external particles. 

Having computed the two-loop amplitude by a genuine perturbative approach 
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without any a priori ansatz on its form, we can investigate whether dual conformal 
properties can be detected even at the level of Feynman diagrams. We have then 
studied the two loop diagrams entering our calculation, out of the mass-shell and in 
three dimensions. Since in three dimensions dual conformal symmetry rules out bub- 
bles, it is immediate to realize that, being some of our diagrams built by bubbles, it 
cannot work at the level of the integrand on every single diagram. A less stringent 
scenario could still allow for the possibility to see dual conformal invariance realized at 
the level of the integrals and after summation of all the contributions. We have made 
many numerical checks but the output is always negative: dual conformal invariance is 
definitively broken at the level of Feynman diagrams. However, this is not in contrast 
with what claimed before, since in three dimensions and in dimensional regularization 
the integrals do not have in general a smooth limit on the mass-shell. 

The organization of the paper is the following. In Section 2 we briefly review the 
ABJ(M) theory in M = 2 superspace formalism. In Section 3 we define light-like 
Wilson loops for the ABJ theory by suitably generalizing the ones for ABJM, and 
compute the two-loop expansion for the four-cusp case. Section 4 contains our main 
results. There, we give details of the one-loop calculation of the four-point scattering 
amplitudes in M > 2 theories and of the two-loop M = 6 ABJ amplitude for four 
chiral superfields. A detailed discussion of our results is given in Section 5 where we 
analyze the amplitude/ WL duality and the related dual conformal invariance of the 
two-loop amplitude, we make a conjecture on its exponentiation and investigate the 
consequences on its strong coupling dual description. Perspectives and open questions 
are also highlighted. Three Appendices fix the notations and provide technical tools 
for carrying on the calculations. 



2. Generalities on ABJ theories 

In three dimensions, we consider M = 2 supersymmetric Chern-Simons theories for 
U(M) x U (N) gauge group, generically coupled to chiral matter. In M = 2 superspace 
the field content is organized into two vector multiplets (V, V) in the adjoint repre- 
sentation of the gauge groups U (M) and U (N) respectively, and four chiral multiplets 
A i and B u i = 1,2 , with A i in the (M, N) and B t in the (M, N) bifundamental 
representations. 

We consider the action (for superspace conventions see Appendix A) 

<S = Sqs + «S mat (2-1) 
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with 



<> cs = J d 3 x d 4 9 jf * dt { ^Tr [VD° {e~ tv D a e tv ) ] + ^Tr [vD° (e~™D a e™ 
5 mat = f d 3 xd 4 9 Tr ( A^tie-v + B^B^ 



+ J d 3 x d 2 9 Tr [h 1 (A 1 B 1 A 2 B 2 ) + h 2 {A 2 B X A X B 2 )] + h.c. (2.2) 

Here hi are two generic complex couplings, while K ± , K 2 are two independent integers 
as required by gauge invariance of the effective action. In the perturbative regime we 
take K 1: K 2 ^> M,N. 

For special values of the h^s we can have enhancement of global symmetries and/or 
R-symmetry with consequent enhancement of supersymmetry [55]. 

Fori^x = — K 2 = K and h\ = —h 2 , the global symmetry becomes U(1)rX SU(2)^x 
SU(2) b and gets enhanced to SU(A) R for [1, 55] 

h x = -h 2 = An/K (2.3) 

For this particular value of the couplings we recover the M = 6 superconformal ABJ 
theory [2] and for N = M the ABJM theory [1]. 

Non-trivial fixed points can be found also for K\ ^ —K 2 [57, 58]. In the large 
M, N limit, choosing 



1 1 



hl = -h 2 = 4tt ^\— 2 + + (2.4) 



we find Af = 2 superCFT's with SU(2) A x SU(2) B global symmetry. These theories 
are dual to string backgrounds where a Romans mass has been turned on [56]. 
More generally, we can take h\ ^ —h 2 satisfying in the large M, N limit 

1^ + 1^ = 32^ + ^ + ^) (2.5) 

This corresponds to a set of M = 2 superCFT's with U(1) A x U(l) B global symmetry. 

The quantization of these theories can be easily carried on in superspace after 
performing gauge fixing (for details, see for instance [57, 58, 59]). In dimensional 
regularization, d = 3 — 2e, and using Landau gauge, this leads to gauge propagators 
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(in configuration and momentum spaces) 



^iD^^i-ft,)^ (2.6) 
A 2 P 



Analogously, the scalar propagators are 



tf a m^r2» - T{1/2 - e) m - 62) 5*5 b 



^'{0,-6,) 5\5 b (2.7) 



where a, b and a, 6 are indices of the (ant i) fundamental representation of U{M) and 
U(N), respectively. 

The vertices needed for loop calculations can be easily read from the action (2.2) 
when expanded in powers of (V, V) up to the desired order. 

The renormalization of the effective action e r = J e s , has been studied up to two 
loops in Refs. [57, 58] using dimensional regularization with dimensional reduction (see 
also [60]). At one loop, there are only finite corrections to the gauge quadratic action 
which in the large M, N limit read 

rSL, = GM " t) / |^M^>^^>) 

f£U = 1] (M - j) / (gi ^M^p^(-p)) (2.8) 

where the G function is defined in (B.l). At two loops, UV divergences appear which 
lead to non-trivial beta-functions. The condition for them to vanish determines the 
conformal fixed points listed in Eqs. (2.3, 2.4, 2.5). 
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3. Light— like Wilson loops for ABJ theories 



The four-cusp Wilson loop for ABJM theory with U (N) x U (N) gauge group proposed 
in [61, 47] is 



{W, 



4/ABJM 



_ 1 fTre i f-/ A,ldz ' i + Tre i ^ A " dz "\ 
2N I J 



(3.1) 



where 7 is a light-like four-polygon closed path. 

The perturbative evaluation has revealed that this expression vanishes at one loop, 
while at two loops it gets non-trivial contributions both from the gauge and the matter 
sectors [47] . In the large N limit one can write 



,(2) 



4/ABJM 



N 2 



{w 4 )P s + (w 4 ; 



matter 



where, up to (9(e) terms, 



(W, 



V(2) 
4/CS 



1 



i ln2 (4^T + (4^T i (fl6 _ 81n2 _ ff2) 



(3.2) 



(3.3) 



as being a constant determined numerically (see Ref. [47]), and 



(W, 



.(2) 

4/ matter 



K 2 



(x 2 3 4Tre^i2 2 ) 2e (x^ 4 4vre^/i 2 ) 2e 



(26) 



+ 



(26) 



--In 5 



1;-! 



X 



21 



7T Z 
T 



(3.4) 



Here 7^ is the Euler constant. Summing the two contributions, one obtains 



(W. 



v(2) 
4/ABJM 



A 5 



( x 13 ^Wl) 2€ 



{2ef 



(^4^H/L) 2 : + I ln 2/^13 ) _ _ ( , 



(26) 



X 



21 



(3.5) 



where A = N/K, ^ 2 WL = 87re 7£ /j 2 and 



C = 3C 2 + 21n2 + 51n 2 2-^ 



(3.6) 



We generalize the definition (3.1) to the larger class of ABJ models, taking into 
account the possibility for the two gauge groups to have different ranks. 

It is convenient to introduce two 't Hooft couplings A = M/K and A = N/K. It 
follows that the perturbative parameter is A = VaT, while 



A- A 



(3.7) 



measures the deviation from the ABJM theory. 
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For gauge group U(M) x U (N) we consider 

W^j = ^mTtumSW + ^Tr u{N) e*^ d *» (3.8) 

Following Ref. [47], it is immediate to realize that this expression is still vanishing at 
one loop, as this result is independent of the choice of the gauge group. 

At two loops, the contributing diagrams are the same as in the ABJM case, but 
with different color coefficients. Generalizing the calculation of [47], in the large M,N 
limit we obtain 

Wil = \(M 2 + iV 2 ) (W A )%> 8 + MN (Walter (3-9) 

(2) (2) 

where (Wi) cs and (yVi) matter are s ^ m given in eqs. (3.3) and (3.4), respectively. 

Inserting their explicit expressions and rescaling the regularization parameter as 

vt* L = #+*' 2 ice"it (3.10) 
up to terms of order e, the final answer reads 

/w\( 2 ) \2 f { x \z^wl) { x 2A^wl) , 1 i 2 ( x \z\ , n i\ 
( ^ 4)abj " M " {2e)^ + 2 ln [Wj > 

(3.11) 

where 

a = fa + 4) < 2 + + 2) ln2 - ?±1 «• + <g2+2) f +10> l" 2 2 (3-12) 
4 8 4 

It is straightforward to check that taking a — > 0, the above expression reduces to the 
result (3.5) for the Wilson loop in ABJM, rescaling (3.10) and constant (3.12) included. 

It is interesting to note that another possible candidate for the Wilson loop in ABJ 
theory is 

<^4>abj = [Tr u{M) e^, A ^ +Tr u{N) e i ^ A ^} (3.13) 

This naturally arises when writing the gauge field as a (M + N) x (M + N) square 
matrix, A = diag(A, A). 

It is easy to realize that this other expression is still one-loop vanishing. At two 
loops, it is given by a slightly different combination of the expressions (3.3, 3.4) 



wSj = I (m3 + N3) m ^ s + { * m2n + n2 m ) W2*t} 

= (M 2 -MN + N 2 ) (W 4 )<§1 + MN (Wetter (3-14) 
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However, provided that we define 

/4 2 L = 2 3+ff2 7re^V (3.15) 
the calculation leads exactly to the same result as before 

( ^ 4)abj = M " — (2^ + 2 ln \4J J 

(3.16) 

where now 

C" = | (a 2 + 2) C 2 + 2(a 2 + 1) In 2 - ^±1 a 6 + (a 2 + l)(a 2 + 5) In 2 2 (3.17) 

Again, taking a — > we are back to the ABJM result (3.5). 

Up to two loops, the two definitions (3.8) and (3.13) for the Wilson loop in ABJ the- 
ory differ only by the choice of the mass scale and the scheme-dependent C constants. 
At this stage we do not have any tool to discriminate between the two. 

4. Four— point scattering amplitudes 

In Chern-Simons matter theories the only non-trivial scattering amplitudes are scalar 
matter amplitudes, as the vector fields are not propagating. In Af = 2 superspace 
language this means having A, B and their complex conjugates as external superfields. 
Given the structure of the vertices, it is straightforward to see that only amplitudes 
with an even number of external legs are non-vanishing. This is consistent with the 
requirement for the amplitudes to be Lorentz and dilatation invariant [36]. 

Each external scalar particle carries an on-shell momentum p a p (p 2 = 0), a SU(2) 
index and color indices corresponding to the two gauge groups. We classify as particles 
the ones carrying (M, N) indices and antiparticles the ones carrying (M, N) indices. 
Therefore, (A 1 , B^) are particles, whereas (B{, Aj) are antiparticles. 

We are interested in the simplest non-trivial amplitudes, that is four-point am- 
plitudes. These are chiral superamplitudes (A l BjA k Bi) and non-chiral superampli- 
tudes (A i A j A k A l ), (BiB j B k B l ), (A i A j B k B l ) plus possible permutations. While for 
the ABJ(M) theories they can all be obtained from (ABAB) by SU(4) R-symmetry 
transformations, for more general Af = 2 models they are independent objects and need 
be computed separately. 

The color indices can be stripped out, as we can write 

A (xi\ Yli z%wli) = E M°<M • • • , *(4)) Ofe) 
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where (X, Z) stay generically for A or B, (Y, W) for B or A and the sum is over 
exchanges of even and odd sites between themselves. We can then restrict to color- 
ordered amplitudes Ai(cr(l), • • • , cr(4)) with a fixed order of the external momenta. 

We compute amplitudes perturbatively, by a direct superspace Feynman diagram 
approach. Precisely, for four-point amplitudes, we evaluate the effective action quartic 
in the scalar matter superfields. Since in a scattering process the external fields are on- 
shell, it is sufficient to evaluate the on-shell effective action. This amounts to require 
the external superfields to satisfy the equations of motion (EOM) 2 

D 2 A i = D 2 Bi = q ? -jf Ai = 13 2 B i = (4.2) 

from which further useful equations follow 

i&^DpA* = id^DpBi = , id^DpAi = id^DpB* = (4.3) 

In principle, setting the external superfields on-shell might cause problems when IR 
divergences appear in loop integrals. We dimensionally regularize these divergences 
working in D = 3 — 2e dimensions, while keeping spinors and tensors strictly in 
three dimensions. We then use the prescription to set the external superfields on-shell 
at finite e. 

To summarize, the general strategy is the following: For a given process and at a 
given order in loops we draw all super-Feynman diagrams with four external scalar su- 
perfields. The corresponding contribution will be the product of a color /combinatorial 
factor times a function of the kinematic variables. We work in the large M, N limit and 
perturbatively in A = MjK\ and A = N/K 2 . To determine the kinematic function, we 
perform D-algebra to reduce superdiagrams to a linear combination of ordinary mo- 
mentum integrals. This is achieved by integrating by parts spinorial derivatives coming 
from vertices and propagators and using the algebra (A. 12) up to the stage where only 
one factor D 2 D for each loop is left. This procedure is highly simplified by the on-shell 
conditions (4.2, 4.3) on the external superfields. We then evaluate momentum integrals 
in dimensional regularization by using standard techniques (Feynman parametrization 
and Mellin-Barnes integrals). 

In momentum space the external superfields carry outgoing momenta (pi,P2, P3, Pa), 
with p\ = and = 0- At the level of the effective action we are allowed to conve- 

niently rename the external momenta, since the p^s are integrated. When evaluating 

2 The actual EOM as derived from the action (2.2) would be D 2 A 1 = -h 1 B 2 A 2 B 1 - h 2 B 1 A 2 B 2 , 
D 2 A 2 = —hiB 1 AxB 2 — h 2 B 2 A\B x plus their hermitian conjugates, and similarly for the B fields. 
However, being us interested in the quartic terms of the effective action, we can safely approximate 
the EOM as in (4.2). 
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the amplitude, the total contribution from every single graph will be given by the sum 
over all possible permutations of the external legs accounting for the different scattering 
channels. 

Mandelstam variables are defined as s = (pi + p 2 ) 2 , t — (pi + P4) 2 , u—(pi + p 3 ) 2 - 
4.1 Scattering at one— loop 

For a generic M = 2 model described by the action (2.2), we first concentrate on the 
chiral amplitudes (A l BjA k B{). 

At tree level and one loop the corresponding contributions are depicted in Fig. 1 
where the four-point interaction comes from the superpotential term in (2.2). 

X 

(a) (6) 

Figure 1: Diagrams contributing to the tree level and one-loop four-point chiral scattering 
amplitude. 

The tree-level amplitudes as coming from Fig. 1(a) are simply given by 

Af ee (A 1 (pi), Bi(p 2 ), A 2 (p 3 ), B 2 (pa)) = h 
A t r(A 1 (p 1 ),B 2 (p 2 ),A 2 (p 3 ),B 1 (p i )) = h 2 (4.4) 

At one loop, we need evaluate diagram 1(6). Performing on-shell D-algebra and going 
to momentum space, the corresponding term in the effective action turns out to be 
proportional to 

/ d 4 6 Tr(A i (jp 1 )B j (p 2 )D a A k (p 3 )D p B l (p 4 )) x /g^ 

i=^i / d 4 6 Tr(A i (p 1 )B j (p 2 )D a A k (p 3 )D l3 B l (p i )) ^'^f (4.5) 

8 J \p 3 + P4r 

where in the second line we have used the results (B.10, B.ll) for the scalar and 
vector-like triangles in dimensional regularization. 

Now, using the on-shellness conditions (4.3), which in the case under exam read 
P3 D a A k (ps) = and p°^DpB\(pi±) = 0, it is easy to see that the final result is zero. 
Since the same pattern occurs for all the permutations of the external momenta, we 
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conclude that the chiral four-point amplitude is one-loop vanishing. This occurs not 
only in the planar limit, but also for any finite value of M, N. 

Notably, the one-loop vanishing of the effective action, quartic in the chiral super- 
fields, can be proved to be true even off-shell [59]. 

We now consider non-chiral amplitudes of the type (A l AjA k Ai). For generic M = 2 
superCFT's we do not expect them to be related to the chiral amplitudes. Therefore, 
a priori there is no reason to expect them to vanish. 

The relevant diagrams for these amplitudes are listed in Fig. 2. 



.4 




Figure 2: One-loop diagrams contributing to non-chiral amplitudes (A l AjA k Ai). 

For each graph we compute the corresponding color /combinatorial factor and per- 
form on-shell D-algebra. We list the results valid for M, N finite (for the time being, 
no large M,N limit is taken). Since we work at the level of the effective action, an 
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overall integral over pi momenta is understood. We also neglect an overall (47r) 2 coming 
from the gauge propagators (2.6). 

Diagram 2(a) : This is the only diagram which involves the chiral interaction vertices 
proportional to h±,h 2 . In this case D-algebra is trivial and the resulting color structure 
gives only double traces. Exploiting the possibility to relabel the integrated momenta, 
the result can be written in a quite compact form 



2(a) = ^ J d 4 9 {(M 2 + \h 2 \ 2 ) Tr(A^p 1 )A l (p 2 ))Tr(A^p i )A J (p 3 )) (4.6) 

+ (hji 2 + h 2 h) TriA'ip^Ajips)) Tr(AJ(p 4 )A t (p 2 )) 

-\h l + h 2 \ 2 Tr(A\p 1 )A^p 2 ))Tr(A\p A )A l (p 3 ))} B( Pl + p 4 ) 



where £>(pi+p 4 ) is the bubble integral defined in (B.8). Repeated indices are understood 
to be summed. 

Diagram 2(6) : In this case the result is a linear combination of single and double 
traces. Single traces are associated with planar graphs and are leading in the large 
M, N limit. D-algebra is easily performed and leads to 



2(6) = J d*6 | - ^ Tr(A*( Pl )Mp 2 )Ai(p,)Mp 3 )) (4.7) 
- TriA^p 1 )A,(p 3 )A^p 4 )A l (p 2 )) 



AK 2 



^ + TF2 ) Tr(A*( Pl )Mp 2 ))Tr(Ai(p 4 )AAP3)) 



Tr^A 1 (pi) Aj(p 3 )) Tr(A j (pi) Ai(p 2 )) } B( Pl + p 2 ) 



Diagram 2(c) : With a convenient choice for the internal momentum, this diagram 
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gives rise to 



2(c) = j d i e\^Tr{D a A\p l )D P A l (p 2 )A^p A )A 3 (p,)) 



N . —8 - 

+ ^ 2 Tr(D a A*(p 1 )A J (p 3 )Ai(p i )D l *A t (p 2 )) 



+ (2^? + 2^f) Tr ( DaAi fa)rf 3 MP2))Tr{A^p A )A 3 {p,)) 

+ Tr(D a A*( Pl )Aj(p 3 )) Tr(Ai(p 4 )D^A t (p 2 )) 

K 1 K 2 



where D-algebra requires integrating two spinorial derivatives on the external fields. 

Using the result (B.ll) for the vector-like triangle in dimensional regularization, 
this contribution vanishes due to the equations of motion (4.3) D a A(pi) = and 
pfDpAfa) = 0. 

Diagram 2(d) : This is the first case where on-shell D-algebra and repeated use of the 
equations of motion allow for a drastic simplification of the final result. We give few 
details of the calculation. 

Computing the color factors we obtain only double trace structures. After per- 
forming D-algebra we are led to 



/ {(2^ + 2^) ^(d-A^p^Am) Tr{pn^)D%^)) 
1 Tr [p^ip^Ajip^j Tr (rP A^(p 4 )D^ 



2 



' (2vr) 3 - 2e k 2 (k + Pl ) 2 (k + Pl + p 4 ) 2 (k-p 2 ) 2 y ' J 



We first integrate by parts the D derivative. Using the equations of motions (4.2,4.3) 
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we obtain 



- Tr (a\ Pi )D & A 3 {p,)) Tr (L>^(p 4 )A(p 2 )) 

(]^ k J)^ A^Q/^y fops 



+ K, 



(27r) 3 " 2e k 2 (k + Pl ) 2 (k + Pl + p 4 ) 2 (k - p 2 ) 2 

]q + 2lq) Tr ( DaAi (PJMP2)) Tr {^{p0 
±_ Tr (^D a A\ Pl )D & A 3 ( P z)) Tr (A (p 4 )A(p 2 )) 

J (2tt)^ k 2 (k + p 1 ) 2 (k + p 1 +p i ) 2 (k-p 2 ) 2 { ' } 

We concentrate on the first integral. The numerator can be rewritten as 

Pi^k^kps = p1 P [k 5l kp a - k 2 C 6a C^] — (k + pi) 2 k 1& - k 2 (k + Pl ) l5 (4.11) 

Now, simplifying the squares at numerator against the ones at denominators we are 
left with a linear combination of scalar and vector-like triangle integrals. Exploiting 
the fact that in dimensional regularization the scalar triangle is zero, while the vector- 
like one is proportional to a bubble integral (see Appendix B), the first term in (4.10) 
reduces to 

/ *° { + Ml ) Tr (^(P0A(P2)) Tr (D^A^p^A^)) (4.12) 

-— Tr [A\ Pl )D 5 A 3 {p z )) Tr (L>^(p 4 ) A(p 2 )) J x ^M^ B( Pl +p 4 ) 

where equations of motion and momentum conservation have been used. Now, inte- 

— 5 

grating by parts the D derivative and using on-shell conditions, it can be further 
simplified to 

/ *° {(^F + 2^) Tr (^(PMM Tr (^(P4)^(ps)) (4-13) 
+ — Tr (A^Ajips)) Tr (^(p 4 )A(p 2 )) J x 8(pi + Pa) 

We can apply the same tricks to the second integral in eq. (4.10). After a bit of algebra, 
we obtain a similar expression which, summed to the rest, leads to the final expression 
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for the box diagram 2(d) in terms of a linear combination of bubbles 

2(d) = / #0 { + ^) Tr (^( Pl )A(P2)) Tr (^(p^-fo)) (4.14) 

+ — ^ Tr (AipJAjfa)) Tr (A>(p 4 )A i (p 2 )) } x [B( Pl + p 2 ) - B( Pl + pA] 

Diagram 2(e) : The result for this diagram reads 

2(e) = - J d 4 9 |^ Tr ^(p^Mp^A^p^A^ (4.15) 

"2^2 Tr (^(Pl)^^(P3)^ Q ^(p 4 )A(p2)) 



Tr (^(pO^A^pa)) Tr (A(p2)T> a ^(p 4 )) 



7 / d 3 ~ 2£ fc k}(k+pi+P3)Sf) 



X (P2)J 



c^ 3 - 2e A; kj_ (k + pi+ p 3 )j 

(27r)3-2e A; 2 (A; + p 1 ) 2 (A; + p 1 +p3) 2 (A;-p 2 ) 2 



Elaborating its numerator, the integral can be rewritten as 

fd^k eC^ + k 1 s (p 1 +p 3 ) S p s 

Jj2^ k^(k + p 1 y(k + p 1+ p 3 nk-p 2 y = C ^(Ps^)HPi+PsU Qv 7 (4.16) 

As proved in Appendix B, the triangle and vector-like box integrals are 0(e) in dimen- 
sional regular izat ion (see eqs. (B.10, B.26)). Therefore, this diagram can be discarded 
when e — > 0. 

Diagram 2(f) : We now consider IP-reducible diagrams. For graph 2(f) using the 
one-loop correction to the gauge propagator given in eq. (2.8) we obtain 

jd A 9 - Tr (^ Q ^(Pi)A(P2)^(p4)^^(p 3 )) (4.17) 

1 ( M- ?p) Tr ^A l (p 2 )D a A^(p 1 )D' 3 A :j (p 3 )A^(p 4 )^ 

Tr (D a A i (p 1 )A i (p 2 )) Tr ^(p 4 %(p 3 )) 



+ K 2 



1 1 2 



AK{ AK 2 KiK- 
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We can integrate by parts the D a derivative. Exploiting the equations of motion, 
the only non-vanishing term is the one where the derivative hits D l3 Aj(p 3 ), giving a 
factor — p?f(p±)ap — ~(P3 + Pi) 2 - By momentum conservation, this cancels against the 
denominator in (4.17) and we finally obtain 



2(/) = - J (n - ?fj Tr {A*( Pl )Mp 2 )Ai(p 4 )AM)) 

+ 7q ( M " t) Tr (MP2)A i (p 1 )Mps)^(p 4 )) 



(4.18) 



+ 



AKf AKi K ± K 2 



Tr (A l (p 1 )A l (p 2 )) Tr (A^p A )A 3 (p^)) \ x B( Pl + p 2 ) 



Diagram 2(g) : Finally, we consider the reducible triangle diagram. Performing on- 
shell D-algebra we produce terms with four spinorial derivatives acting on the external 
fields. After integrating by parts one of these derivatives, using on-shell conditions and 
momentum conservation and relabeling internal and external momenta, we can write 
the result as 



-J d A 9 ly-j^Tr (a^D^AMD-A^pM^ 
~J^K 2 Tr ( Ai (pMj(P3)D a A^(p 4 )D A i (p 2 )) 



+ 



x 



+ Tr (A\ Pl )D P A^)) Tr(L>^(p 4 )^(p 3 )) 



(pi+p 2 ) 7 a f d 3 ~ 2e k (k + Pi y s k 5 p 



(P1+P2) 2 7(2tt) 3 -^ k 2 (k+ Pl ) 2 (k-p 2 ) 2 
We can elaborate the numerator of the integrand to obtain 



d^k k 2 5l+pfksp 



(2tt) 3 - 2£ k 2 (k + Pl ) 2 (k-p 2 y 



SlB( Pl +p 2 )+ P J s % 



V8/3 



B( Pl +p 2 ) 



51- 



PT(.P2)5I3 
(Pi +P2Y 



(4.19) 



(4.20) 



where eq. (B.ll) has been used together with on-shell conditions. Inserting back into 
eq. (4.19), observing that on-shell (pi + p 2 ) ya pJ S (p 2 )s/3 = (Pi + P2) 2 (P2)a/3 where (p 2 ) a p 
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vanishes when contracted with D Ai(p 2 ), we obtain 



Tr (A i {p 1 )I^A i {p 2 )D a AHp^Aj{p z ] 



iK 2 
M 



+ 



\s + J£l) Tr (A\ Pl )D P A^)) Tr{D«A\p,)A 3 (p,)) 



K 



(Pl)a/3 



x B(p 1+ p 2 ) (4.21) 

(Pi +P2) 2 



On-shell conditions are once again helpful for reducing the structure of spinorial deriva- 
tives acting on the external fields. In fact, integrating by parts the D a derivative we 
produce a term p^ that, contracted with (pi) a p, cancels (pi +p 2 ) 2 at the denominator. 
We finally obtain 



!(<?) = - / d 4 j-^- Tr {A^AM^ipMjips)) (4.22) 



2 



+ {lq + ^f) Tr ( Ai (P^MP2)) Tr (AHp^Ajfaj) \ x B( Pl +p 2 ) 



We are now ready to sum all the results and obtain the one-loop effective action 
needed for the evaluation of (A l AjA k Ai) amplitudes. 

Having reduced all the expressions to strings of external superfields with no deriva- 
tives acting on them, multiplied by bubble integrals, we can group them according to 
their trace structure. We have single trace contributions from diagrams 2(6), (/), (g) 
and double trace contributions from 2(a), (6), (d), (/), (g). Collecting them all, we ob- 
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tain 



J d 4 9 | - Tr{A\ Pl )MP2)A^p,)A 3 {p,)) A (±- + ±^ x + p 2 ) 
- rr(^( Pl )^(P3)^(p4)^(p 2 )) A (±- + ^ x B(p! + P2 ) 



(4.23) 



+ Tr(^( Pl )A(p 2 )) Tr(^(p4)^(p 3 )) * 



1 f\hl\ 2 + \h^ 2 



2 V 16tt 2 



+ Tr{A\ Pl )A 3 { Pz )) Tr(^( P4 )A(p 2 )) - x Bfa + p 4 ) 

Tr(^(p 1 )A l (p 2 ))Tr(A l (p 4 )^(p3)) + x B( Pl + p A ) 



327T 2 

First of all, we observe that for M, N finite the quartic effective action, and consequently 
the amplitude, vanishes when K 2 = —K 1 and h 2 = —hi = Att/K. This is exactly the 
M = 6 superconformal fixed point corresponding to the ABJ theory. This result was 
expected and provides a non-trivial check of our calculation. In fact, in the ABJ model 
the non-chiral amplitude is related to the chiral one by SU (4) symmetry and we have 
already checked that the chiral amplitude is one-loop vanishing. 

Taking M, N large and assuming the hi couplings of order of 1/K i: only single trace 
contributions survive in (4.23). In this case, the amplitude will vanish for K 2 = —Ki, 
independently of the values of the chiral couplings. In particular, we have a vanishing 
non-chiral amplitude for the whole set of M = 2 superCFT's given by the condition 
(see eq. (2.5)) 

M 2 + N 2 = ^ (4.24) 

We observe that the amplitude never vanishes for theories with K 2 ^ —Ki, in particular 
for superCFT's which correspond to turning on a Romans mass in the dual supergravity 
background. 

The same pattern occurs for the [Bi& B^B 1 ) amplitudes. In fact, repeating the 
previous calculation we obtain exactly the same expression (4.23) as a consequence of 
the Z 2 symmetry of the action under exchange V -H- V, A -h- B, M -h- N, Ki -h- K 2 
and hi ■<-)■ h 2 . 

Finally, we need consider mixed amplitudes of the type (A l AjB k Bi). The con- 
tributing diagrams are still the ones drawn in Fig. (2) with obvious substitution of one 
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(A, A) couple with a (B, B) couple. Applying the same procedure as before, we obtain 



For M,N finite, these amplitudes vanish only at the ABJ(M) fixed point (2.3), as 
expected. However, in contrast with the previous case, in the large M, N limit a non- 
trivial dependence on the chiral couplings survives, which restricts the set of superCFT's 
with vanishing one-loop amplitudes to be only the ABJ(M) models. 

In the ABJM case, this result is consistent with what has been found in Ref. [46] 
in components and massive regularization, and in [51] by means of the generalized 
unitarity cuts method. 

4.2 Scattering at two loops 

We restrict to the ABJ model (2.3) for which we have found that all the four-point 
amplitudes vanish at one loop. This result is consistent with the one-loop vanishing 
of the Wilson loop and leads to conjecture the existence of a WL/amplitudes duality 
for this theory. To give evidence to this conjecture, we evaluate four-point scattering 
amplitudes at two loops. 

We study amplitudes of the type (A l BjA k Bi), where the external A, B particles 
carry outgoing momenta p±, . . . ,p± (pf = 0). 




(4.25) 



+ 2Tr(A i 



+ 2Tr(A i 




x B{pi + p 4 ) 
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At two loops, in the planar sector, the amplitude can be read from the single trace 
part of the two-loop effective superpotential 



T^[A,B] = J d 2 6d i Pl ...d i PA (2nf5^\Y J Pi) * 

^e^'tr (^( Pl )^(p 2 )A fc (p3)^(p 4 )) J2 M^ x \ Pl , . . . ,p 4 ) (4.26) 

X=a 

where the sum runs over the six 1P1 diagrams in Fig. 3, plus the contribution from 
the IP-reducible (1PR) graph in Fig. 3(g) where the bubble indicates the two-loop 
correction to the chiral propagator. 

In (4.26) we have factorized the tree level expression, so that M {x) (pi,... ,P4) are 
contributions to A^/A*™ 6 . 

In order to evaluate the diagrams we fix the convention for the upper-left leg 
to carry outgoing momentum pi and name the other legs counterclockwise. The 
momentum-dependent contributions in (4.26) are the product of a combinatorial fac- 
tor times a sum of ordinary Feynman momentum integrals arising after performing 
D-algebra on each supergraph (more details can be found in [59]). There are a total of 
four diagrams of the classes (6), (c), (/) and (g), eight diagrams of the classes (d) and 
(e) and two diagrams of the class (a). The color/flavor factors C; for them are given by 



Ca = (47r) 2 ^ Cb = (47r) 2 ^ C c = (4tt) 2 *£=£=k 
C d = (4tt) 2 C e = (4tt) 2 AA C f = -(4vr) 2 AA (4.27) 

while diagram (g) contains subdiagrams with different color/flavor factors and these 
cannot be factorized. 

Diagram 3(a) : We begin with the simplest graph which, after D-algebra, reduces to 
the following factorized Feynman integral 

V s = *e f d^J^ -{P1+P2? , . 

a * J ( 27r )d(2ir) d k 2 (k+ Pl +p 2 yP(l-p 3 -p 4 ) 2 { - } 

where \x is the mass scale of dimensional regularization. 

The k and the I bubble integrals can be separately evaluated using the result (B.8), 
so obtaining 

K = ~G[1, l] 2 ( ^ (4.29) 



2 



x s r 

In order to obtain the corresponding contribution to the amplitude, we need sum over all 
the independent configurations of the external momenta. Inserting the corresponding 



- 22 - 






(a) 







(d) 



Figure 3: Diagrams contributing to the two-loop four-point scattering amplitude. The 
dark-gray blob represents one-loop corrections and the light-gray blob two-loop ones. 



color/flavor factors we obtain 



M {a) = -8tt 2 (A 2 + A 2 )G[1,1] S 



2\ 2e 



+ 



2\ 2e N 



-C 2 (A 2 + A 2 ) + 0(e) (4.30) 



Diagram 3(6) : After D-algebra, it gives 

d d k d d l 



2(p 3 + Pif 



(27r) d (27r) d P (I + k) 2 (k - p 4 ) 2 (k + p 3 f 



(4.31) 



Performing the I integral with the help of Eq. (B.8) , we obtain a triangle integral with 
a modified exponent in one of its propagators 



d d k 



2(p 3 + p 4 )' 



(2vr) d {k 2 y/ 2 +t (k - Pa) 2 (k + p 3 f 



(4.32) 



We Feynman-parametrize the denominator and integrate over the momentum k. Tak- 
ing into account that we are working on-shell (pf = 0) we finally obtain 



up 



/i 4e 2sG[l, 1] T(l + 2e) f dMMPs A - 



1/2+e 



(47r) d / 2 r(l/2 + e) J (f3 1 (3 2 pj + (3 1 (3 3 p 2 + f3 2 (3 3 s) 1+2e 
o 

pL->o 2G[l,l]r(l + 2e)r 2 (-2e) //i 2N 2 ' 



(47r) d / 2 r(l/2 - 3e) V s 
where T 2 (— 2e) signals the presence of an on-shell IR divergence. 



(4.33) 
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Summing over the four inequivalent configurations of the external legs multiplied 
by the correct vertex factors, the contribution to the amplitude reads 



M {b) = 8tt 2 (A 2 + A 2 ) 



A 2 + A 2 



(2c) s 



G[l,l]r(l + 2e)r 2 (-2e) 
(4vr) d / 2 r(l/2 - 3e) 

+ 



7re 



(2eY 




2c x 



2(2 + 0(6) 



(4.34) 



Diagram 3(c) : This diagram may result problematic, being infrared divergent even 
when evaluated off-shell. In fact, after D-algebra, the particular diagram drawn in 
Fig. 3(c) gives rise to the following integral 

= f d d k d d l Tr{ lt g vlpl(T )^{k-p z Y {k + p^Ypl 

2 J {2ir) d {2n) d I 2 (I + k) 2 k 2 (k - p A ) 2 (k + p 3 ) 2 1 ' ' 

which, performing the bubble I integral and using the identity 

Tr{ llxlvlpla )rt{k-pzY '(k+PiYPl = (P3+PA) 2 k 2 -pl(k + p A ) 2 -pl(k-p 3 ) 2 (4.36) 
can be separated into three pieces 

K 1 = \ X?f - \C[1, l]G[l, 3/2 + e] {pi)' 2 ' - \g[1, 1]G[1, 3/2 + e] (p 2 )" 2e (4.37) 

While the first term is the off-shell well-behaving Feynman integral in Eq. (4.33) that 
in the on-shell limit produces 1/e poles, the second and third terms are badly divergent 
even off-shell. However, we can show that these unphysical divergences are cured when 
we add the 1PR diagrams corresponding to two-loop self-energy corrections to the 
superpotential, as depicted in Fig. 3(g). 

In fact, for example the contribution from the diagram with the two-loop correction 
on the fourth leg as drawn in the picture, yields 

V 4 g = -8tt 2 (8AA - A 2 - A 2 ) G[l, 1] G[l, 3/2 + e] (p 2 4 )- 2t + 32tt 2 AA p\ B(p 4 ) 2 (4.38) 

where color factors have been included. 

It is easy to realize that the first term of this expression is off-shell infrared di- 
vergent, but precisely cancels the third term in (4.37) when all the vertex factors of 
diagram 3(c) are taken into account. On the other hand, the second term in (4.38) 
comes from a double factorized bubble which vanishes on-shell. 
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Since in a similar way the second term in (4.37) gets canceled by a diagram with a 
two-loop correction on the third leg, summing diagrams 3(c), (g) and their permutations 
we are finally led to the following interesting identity 



(4.39) 



A 2 + A 2 

Diagram 3(d) : Diagrams of type (d) may be calculated using Mellin-Barnes techniques. 
Specifically, after D-algebra the diagram in figure gives rise to the integral 

d d k d d l Tr(^ l/ ^ a )p1(p 3 +p 4 Y(k+p 4 y(l-p 4 y 



(2ix) d (2vr) d (k + p 4 ) 2 (k - p 3 ) 2 (k + I) 2 (I - p 4 ) 2 1 2 
Using the identity (B.6) and the on-shell conditions, it can be rewritten as 



(4.40) 



-sT(l/2-e) j" dz 

(47r) d / 2 r(l - 2e) J 2t7i 

—ioo 

d d k 



X fi 



Ac 



r(-z|l + z|3/2 + e + z| - 1/2 - e - z) 
1 



(4.41) 



(2ll) d (P)3/2+,+z p + p4 )2]-z (fc _ p3 )2 

The integral over can be easily performed by Feynman parametrization, leading to 
4e r(l/2-e|l+2e| -2e) 



2tf = 



(47r) d r(l-2e|l/2-3e) 

ioo 

/ ^r(l + z|3/2 + e + z| - 1/2- 

—too 

r 3 (l/2-e)r(l + 2e)r 2 (-2e) 
(47r) d r 2 (l - 2e)r(l/2 - 3e) (s//x 2 ) 2e 



2e 



(4.42) 



where the remaining integral over the complex variable z has been performed by using 
the Barnes first Lemma (B.7). 

Taking into account the eight permutations with corresponding flavor/color factors 
we obtain 



M®= -16tt 2 (A 2 + A 2 ) 



r 3 (l/2-e)r(l + 2e)r 2 (-2e) 
(47r) d r 2 (l - 2e)r(l/2 - 3e) 



A 2 + A 2 




( s ^ 




1 1 ) 


4 


(2e) 2 


y47re-^ E fi 2 J 


(2e) 2 





-2e 



2 C 2 + 0(e) 



(4.43) 
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Diagram 3(e) : Using the identities derived in [59] it is possible to write this diagram 
as a combination of diagrams (b) and (d). It holds that 



(1 + S^Vf = (1 + S z4 )Vf + Vf - p\B{p z ) 2 - p 2 A B{p,) 2 (4.44) 

where T> s e l is the particular diagram of type (e) drawn in the figure and the operator 
(1 + 1S34) symmetrizes the diagram with respect to the third and fourth leg. Notice the 
presence of a double factorized bubble which can be dropped on-shell. Accounting for 
all permutations and flavor/color factors we thus find that 

M {e) = 4 2 A \ 2 (M (d) + 2M (h) ) (4.45) 
A^ -|- A^ 

Diagram 3(/) : The most complicated contribution comes from this diagram, as it 
involves a non-trivial function of the s/t ratio. Surprisingly, after some cancelations it 
turns out to be finite. 

The D-algebra for the specific choice of external momenta as in figure results in 
the Feynman integral 

^234 = ;y 4e f d d k d d l -Tr{ lillvlpla )p^ 2 kn° 

f ^ J (2ir) d (2Tr) d ki(k-p 2 y(k + i + p 3 y(i-p 4 yp { - 1 

Again, using Eq. (B.6) for the k integral and working on-shell, we obtain 



p234 = r(l/2 - e) 



f " (4^)^(1 - 2e) 



ICO 

[ ^r(-z)r(i + z )r(3/2 + e + z)r(-i/2 - e - z) 

J 2m 



-too 



x , Ae [ d d l %Wp7j P>2 (* + PsY 1° Um 

^ J (27T) d P (I - p,Y [(/ + P3 ) 2 ]" Z [(I +P2+ P3 ) 2 ] 3/2+e + Z ^ ' 1 

It is convenient to separate the / integral in two pieces by using the on-shell identity 

Tri^^p^p^l+p.YF = -(s + t)l 2 + Tr{ 1 , lvlpla )p l ip^r (4.48) 

on-shell 

The first piece in (4.48) contains an I 2 factor which cancels the I 2 propagator in (4.47) 
leading to a simple triangle which is straightforwardly evaluated as we did for the 
previous diagram. 

The second piece is a vector-box integral which after Feynman parametrization 
can be written in terms of a second 1-fold Mellin-Barnes integral. Interchanging the 
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order of the two Mellin-Barnes integrals and solving for the original one in (4.47) with 
the first Barnes lemma (B.7), the total result is 



Vf 4 



( s + t)r 3 (i/2-6)/^ 

(47r) d r(l/2-3e)r 2 (l-2e) 



r(l + 2e)r 2 (-2e) 



; l+2e 



+ 



loo 

+ ^/ ^ r (- v ) r (- 2e - v ) r (- 1 - 2e - v ) r2 ( 1 + v ) r (2 + 2e+v) (^) V 



(4.49) 



The contour of the Mellin-Barnes integral in the second term of the last expression is 
not well-defined in the limit e — > 0, reflecting the presence of poles in e. The reason is 
that in this limit the first pole of T(— 1 — 2e— v) collapses with the first pole of T 2 (l +v). 
In order to have a well defined contour in the e — > limit, we can deform the contour 
so that it passes on the right of the point v = — 1 — 2e and include the residue of the 
integrand in this point. Surprisingly, it turns out that this residue exactly cancels the 
first term in (4.49) so that we obtain a simple 1-fold Mellin-Barnes integral which is 
finite in the limit e — > 



Vf 4 



(i + s /t)r 3 (i/2- e ) 



(4.50) 



+ioo 



X 



(47r)<T 2 (l - 2e)r(l/2 - 3e)(t/^ 2 ) 2e 



This integral can be calculated in the e — > limit by closing the contour and performing 
the infinite sum of all the residues inside it. Taking into account all four permutations 
of the diagram and flavor/color factors we finally obtain 



M (f) = \\ (| ln 2 (s/t) + 3C 2 ) + 0(e) 



(4.52) 



We are now ready to collect the partial results (4.30, 4.34, 4.39, 4.43, 4.45, 4.52) 
and find the four-point chiral superamplitude at two loops. After some algebra, and 
redefining the mass scale as 



the result can be cast into the following compact form 




(4.53) 



(4.54) 
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where A = vAT = \JMN / k and C4 is a constant given by 

C A = (4 - fa 2 ) C 2 + (1 + \o 2 ) (3 + ±a 2 ) In 2 2 (4.55) 

We note that the mass scale and the constant depend non-trivially on the parity- 
violating parameter a defined in eq. (3.7). Since it only appears as a square, parity is 
not violated at this stage. As a check, we observe that for a = the result reduces to 
the ABJM amplitudes computed in [52]. 

5. Discussion 

We now discuss the main properties of our result (4.54) for the four-point amplitude 
at two loops. 

First of all, in the ABJM case (er = 0) the result coincides with the one in [51] 
obtained by applying generalized unitarity methods. In particular, the effective mass 
scale is the same and the analytical expression for the constant C A \ a ^. = 4£ 2 + 3 In 2 2 
matches the numerical result of [51]. 

In Ref. [51] the result has been found by assuming a priori that in the planar limit 
dual conformal invariance should work also at quantum level. In fact, an ansatz has 
been made on the general structure of the amplitude which turns out to be a linear 
combination of integrals that, if extended off-shell, are well defined in three dimensions 
and exhibit conformal invariance in the dual x-space (pi = Xi — Xj+i). 

On the other hand, our calculation relies on a standard Feynman diagram approach 
which does not make use of any assumption. The identification of the two results is 
then a remarkable proof of the validity of on-shell dual conformal invariance for this 
kind of theories. 

5.1 Amplitudes/ WL duality 

In the general ABJ case, if we write the Mandelstam variables in terms of the dual 
ones, s = rr 2 3 and t = x\ A , up to a (scheme-dependent) constant our result matches 
those in Eqs. (3.11, 3.16) for the two-loop expansion of a light-like Wilson loop, once 
we formally identify the IR and UV rescaled regulators of the scattering amplitude and 
the Wilson loops as fi' 2 = I/^wl or f 1 ' 2 — ^-IflwL- 

Since the Wilson loop is conformally invariant in the ordinary configuration space, 
the identification of the two-loop amplitude with the corresponding term in the WL 
expansion is a further proof of dual conformal invariance in the on-shell sector of the 
theory. 

We remind that the two results (3.11, 3.16) in Section 3 correspond to two possible 
definitions of Wilson loop in ABJ models. At this stage, the result (4.54) seems to 
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match both. However, if we compare the rescaled mass regulators, we see that apart 
from the sign of the Euler constant, in the result for the amplitude // looks like \j! wl 
in eq. (3.10), whereas it is quite different from ii'{ VL in Eq. (3.15). Although there 
is no particular reason for the two mass scales to match exactly, this might be a first 
indication that the definition (3.8) for the light-like Wilson loop dual to scattering 
amplitudes is preferable. 



5.2 Dual conformal invariance 

As for the M = 4 SYM case, in the ABJ models the two-loop on-shell amplitude 
divided by its tree-level contribution, when written in terms of dual variables has the 
same functional structure as the second order expansion of a light-like Wilson loop. 
Wilson loops are invariant under the transformations of the standard conformal group of 
the ABJ theory, even though UV divergences break this symmetry anomalously. Hence, 
the on-shell amplitude should inherit this symmetry, possibly anomalously broken by 
IR divergences. 

In fact, in the M = 4 SYM case where the amplitudes/ WL duality also works, the 
perturbative results for planar MHV scattering amplitudes divided by their tree-level 
contribution can be expressed as linear combinations of scalar integrals that are off- 
shell finite and dual conformally invariant [62, 5] in four dimensions. Precisely, once 
written in terms of dual variables, the integrands times the measure are invariant under 
translations, rotations, dilatations and special conformal transformations. Dual con- 
formal invariance is broken on-shell by IR divergences that require introducing a mass 
regulator. Therefore, conformal Ward identities acquire an anomalous contribution [8] . 

A natural consequence of our findings is that the two-loop result for three di- 
mensional ABJ(M) models should also exhibit dual conformal invariance, and then it 
should be possible to rewrite the final expression (4.54) for the on-shell amplitude as a 
linear combination of scalar integrals which are off-shell finite in three dimensions and 
manifestly dual conformally invariant at the level of the integrands. Indeed, for the 
ABJM case this has been proved in Ref. [51] where the amplitude has been obtained 
by unitarity cuts method, based on the ansatz for the amplitude to be dual conformal 
invariant. 

Following [51], we introduce a set of independent scalar integrals h s , hs, hs, hs = 
hs—hs which correspond to the following off-shell, three dimensional, dual conformally 
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invariant expressions 
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plus their t-counterparts obtained by cyclic permutation of the (1, 2, 3, 4) indices. Their 
graphical representation is given in Fig. 4. 
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Figure 4: Graphical representation of dual conformally invariant integrals. 

The appearance of the particular combination hs — J 4s is not an accident. In fact, 
due to the presence of internal cubic vertices, the integrals hs, h s are IR divergent also 
off-shell and then ill-defined in three dimensions. Dual conformal invariance would 
require to discharge these integrals. However, as we show in Appendix C, taking the 
linear combination hs — 7 4s the off-shell divergences cancel and h s is well-defined in 
three dimensions. 

The on-shell evaluation of these integrals in D = 3 — 2e dimensions reveals that 



'2s 



0{e 2 ) 



'3s 



'3i 



-hs ~ I 



It 



(5.5) 
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Thus, the actual basis for the four-point scattering amplitude reduces to h s , I$ s , h t , ht- 
Evaluating them and defining Jj 2 = 87re~ 7E /i 2 , one finds [51] 3 



Iu + ht 



1 



hs + h 



5t 



167T 2 
1 

"8^2 



(s/fi 2 ) 



2e 



+ 



+ 



m 2 r 

2e 

{t/n 2 )' 2 ' 



+ 2 — 21n2 

{sin 2 )' 2 ' 



(2e) 2 (2e) 2 2e 

(0 +21n2-2-31n 2 2-4( 2 



-iln 

2 



+ 0(e) 

m 2 r 2 ' 

2e 

0(e) 



(5.6) 



(5.7) 



Using these results, it is easy to see that for the ABJM theory, the two-loop amplitude 
(4.54) for er = can be written as 



Mf 



ABJM 



(4ttA) 2 



-h s + I l8 + (S <H> t) 



(5.8) 



Remarkably, this linear combination not only reproduces correctly the non-trivial part 
of the amplitude, but also fits the numerical constant Cj±. In particular, it is such that 
the non-maximal transcendentality terms in (5.6), (5.7) cancel. 

We can now generalize this analysis to the ABJ models where the amplitude has 
the same functional structure of the ABJM one, except for a non-trivial dependence on 
the a parameter in the mass scale and in the constant. Because of the appearance 
of <7, we find that in terms of the integrals (5.6, 5.7) given as functions of the Jx 2 scale, 
the M. 4 ratio can be written as 

2 



MP 



ABJ 



a 



-hs(fi ) + 1 + In 2 I ls (fi 2 ) + (s<*t) 



where C rcs is a residual constant given by 



C res = A 2 cr 2 I In 2 2- ^2 + hi2 



(5.9) 



(5.10) 



The non-trivial appearance of a 2 in the coefficients might reflect the fact that in the 
ABJ case one cannot factorize completely the color dependence outside the combination 
of integrals. This would suggest that a generalization of the unitarity cuts method 
should be employed where the trace structures are not stripped out. Nevertheless, it 
is not difficult to see that the application of such a method would never reproduce the 
In 2 coefficient in front of hsiht- 



3 In [51], the constant part of (/5 S + Iu) has been evaluated numerically. However, a posteriori one 
can check that the numerical factor is well reproduced by the analytical expression in (5.7). 
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One would be tempted to conclude that ABJ amplitudes cannot be computed by 
unitarity cuts methods. However, a way out is to start from a linear combination of 
(5.6, 5.7) integrals where the mass parameter has been rescaled as fi 2 — > fi' 2 = p 2 2' j2 / 2 . 
Doing that, we find that the two loop ratio can now be written as 



ABJ 



= (4ttA) : 



-J 5s (// 2 ) + / ls (/i /2 ) + ( S ^t) 



+ CL 



where 



C' = \ 2 a 2 



2 + 



a 



In 2 2 



<2 



(5.11) 



(5.12) 



4 J "' 4 " 

The situation has drastically improved, since rational coefficients in front of the integrals 
indicate that the same result could be obtained by unitarity cuts method. However, in 
that approach the question of why and how fixing a priori a non-standard mass scale 
in the dual invariant integrals remains an open problem. 

Except for the particular case M = N, in general the basis of scalar integrals 
selected by dual conformal symmetry reproduces the four-point amplitude only up to 
a constant. This is a quite different result compared to what happens in the ABJM 
and M = 4 SYM cases where dual conformal integrals reproduce exactly the four-point 
amplitude. However, at the order we are working, the difference is only by a constant 
and dual conformal invariance is safe, as well as the anomalous Ward identities which 
follow. 

At higher loops, we expect the non-trivial dependence on a to affect also the terms 
depending on the kinematic variables. It would be very interesting to check whether 
this phenomenon may spoil dual conformal invariance or higher order amplitudes could 
still be expressed as (ex-dependent) combinations of dual conformal invariant integrals. 
For this reason, it would extremely important to evaluate the amplitude at four loops. 

In the ABJM theories, comparing the result for the amplitude obtained by ordinary 
perturbative methods with the one obtained by unitarity cuts method, we can write 



D = 3 - 2e 
on — shell 



D = 3 - 2e 
on — shell 



(5.13) 



where Jj are momentum integrals associated to the Feynman diagrams in Fig. 3, 
whereas /„ are the scalar integrals (5.1-5.4). 

Since the linear combination on the right hand side, when written strictly in D = 3 
with x 2 i+1 7^ is dual conformally invariant, the natural question which arises is 
whether also the left hand side shares the same property. 

In order to answer this question, we investigate the behavior of the Feynman inte- 
grals Ji under dual conformal transformations, off-shell and in three dimensions. After 
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,2 



rewriting them in terms of dual space variables, we implement the inversion — > -^i 

x i x j 

and d d Xi — > f-rh, which is the only non-trivial conformal transformation to be checked. 



As in four dimensions the invariance under inversion rules out triangle and bubble- 
like diagrams, similarly in three dimensions it forbids the appearance of bubbles. There- 
fore, just looking at the integrands, we see that the integrals associated to diagrams 
3(a)-3(o) cannot be separately dual conformal invariant. Moreover, despite the fact 
that diagrams 3(d)-3(/) consist of triangles only, non-trivial numerators spoil invari- 
ance under inversion, as well. 

Nevertheless, these considerations on the integrands may fail in very special cases. 
As an example, we consider the double bubble diagram 3(a). In dual coordinates, the 
corresponding integral reads 



Performing inversion, this integral gets mapped into a double triangle integral 



If we evaluate them off-shell, we obtain B = T = 1/64 (see eqs. (B.8, B.9)) . Therefore, 
the double bubble diagram is invariant under inversion at the level of the integral, even 
if it is not invariant at the level of the integrand. 

Motivated by this example we may wonder whether dual conformal invariance on 
the left hand side of eq. (5.13) could be restored at the level of the integrals. By 
numerical evaluating the integrals associated to the independent topologies 3(6), 3(d) 
and 3(/) and to their duals, obtained by acting with conformal inversion, we find that 
every single integral is not by itself dual conformally invariant. 

However, one may still doubt that the situation could improve when summing 
over all scattering channels, or summing all the contributions to get the total off-shell 
amplitude. 

We find that, even if for every single diagram the sum over permutations of external 
legs definitively improves the result, as for a large sample of momentum configurations 
the integral and its dual look very close to each other, they never appear to be exactly 
invariant, neither does the total sum. 

We thus conclude that the off-shell amplitude computed by Feynman diagrams is 
not dual conformally invariant. In other words, the identity (5.13) is not an algebraic 
relation between different basis of integrals, as if it were the case it should be valid 
for any value of the kinematic variables. Instead, it holds only when the integrals are 




(5.14) 




(5.15) 
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evaluated on the mass-shell and in dimensional regularization. This is not puzzling if 
we take into account that in three dimensions and in dimensional regularization the on- 
shell limit is not a smooth limit for the integrals. It would be interesting to investigate 
what happens when using a different regularization, for example the one suggested in 
[63, 64]. 

In any case, our result reinforces the statement that dual conformal invariance is a 
(anomalous) symmetry only of the on-shell sector of the theory. 

5.3 BDS-like ansatz 

The striking correspondence between the four-point amplitudes of ABJM theory at 
two loops and the one of 4d M = 4 SYM at one loop led us to conjecture [52] (see also 
[51]) that a BDS-like ansatz [53] may be formulated also for the three-dimensional case 

Mi = e ^+^(ln 2 (f)+8C 2 +6 1n2 2 )+C(A) ^ ^ 

where fcs(X) is the scaling function of ABJM. The analogy between the two theories 
is due to the fact that they share similar integrable structures, with asymptotic Bethe 
equations related by an unknown function h(X) [19]- [22]. This leads to a connection be- 
tween anomalous dimensions of composite operators and, in particular, to the following 
relation between the scaling functions [19] 



fcs(X) = ^= 4 (A) 



(5.17) 



in terms of the interpolating function h(X) that needs to be determined in perturbation 
theory. Since the scaling function governs the coefficients of the kinematic part of the 
four-point amplitude in M = 4 SYM by means of the BDS exponentiation, one is 
tempted to conjecture that an analogue resummation may also hold in the ABJM case, 
giving rise to equation (5.16). This formula is confirmed at two loops by the results of 
[52, 51]. Since at weak coupling h(X) is known up to the forth order [31, 32, 33], we 
easily find 

fcs(X) = AX 2 - 24 C 2 A 4 + 0(X 6 ) (5.18) 

and the ansatz (5.16) provides a prediction for the four-loop expression of the finite 
remainder (in the notation of [53]) for the ABJM four-point scattering amplitude 
[52] 

F 4 < 4 > = £ In 4 Q + A 4 Q In 2 2 - In 2 (0 + Consts (5.19) 

Now we discuss how this scenario might be affected by the generalization to the 
ABJ case, where integrability is not expected to be trivially preserved. We first note 
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that the final expression (4.54) is the result of summing many contributions which in 
general are proportional to the homogeneous couplings A 2 , A and to the mixed A A 
one. It is interesting to observe that at this order it is possible to redefine jj? in such 
a way that in all the terms depending on the kinematic variables, the contributions 
proportional to the homogeneous couplings cancel, leading to an expression which is 
basically identical to the one for ABJM, except for the substitution N 2 — y MN. This 
also happens for the WL computed in Section 2. 

This phenomenon has been also observed in the evaluation of the spin-chain Hamil- 
tonian associated to the two-loop anomalous dimension matrix for single-trace opera- 
tors [65] and in the two-loop contribution to the dispersion relation of magnons [18]. 
This special dependence on the coupling constants is a signal that parity symmetry 
along with integrability are preserved at least at two-loop order even if ABJ theory is 
manifestly parity breaking. At four loops only the dispersion relation, i.e. the eigen- 
value of the Hamiltonian of spin chains with a single excitation, is known to date. 
It has been confirmed by explicit computations and the following expansion for the 
interpolating function has been found [31, 32, 33] 

h 2 (\, a) = \ 2 - A 4 [4C 2 + C2 (J 2 ] (5.20) 

From the last term of this expression it is clear that, even if the departure from the 
ABJM case becomes non-trivial, still the function turns out to depend quadratically 
on a and thus parity breaking is not visible. This might indicate that integrability is 
not broken also at the four-loop level. 

Therefore it seems plausible that, at least up to four loops, the planar limit of 
ABJ scattering amplitude could behave in the same way as in the ABJM case, being 
governed by the scaling function /(A, a) obtained through (5.17) where we replace h(\) 
with h{\, a). 

Since at weak coupling h{\, a) is known up to four loops [31, 32, 33], we find 

f cs (\ <t) = 4A 2 - 4 (6 + a 2 ) (2 A 4 + 0{\ & ) (5.21) 

It is then easy to see that the four-point amplitude at order A 2 can be identified with the 
first order expansion of an exponential of the type in eq. (5.16) with /cs(A, cr) = 4A 2 . 

Moreover, this would lead to a prediction for the four-loop expression for the finite 
remainder of the ABJ four-point amplitude to be given by 

= ^ In 4 (0 + A 4 [1 (1 + \a 2 ) (3 + \a 2 ) In 2 2 - (l + ^ 2 ) (2] ^ (|) + Consts 

(5.22) 
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At four loops the parity-violating parameter is expected to play an active role and 
the theory could present a very different behaviour with respect to the ABJM case. It 
would be interesting to check this expression by a direct computation. 

5.4 The amplitude at strong coupling 

In M = 4 SYM a recipe for computing scattering amplitudes at strong coupling has 
been proposed by Alday and Maldacena [38] within the context of the AdS/CFT cor- 
respondence. According to their prescription, the amplitude for n gluons is obtained 
by computing the minimal area of a surface in the AdSs dual background, ending on a 
light-like n-polygon, whose edges are determined by the gluon momenta 

R AdS* . / 1 \ 

M = e~^^ A + O l-j=) (5.23) 

Here RAds 5 stems for the AdS radius and determines the dependence on the 't Hooft 
coupling (R 2 AdS5 = \fg 2 N = V^sym)- 

After a suitable regularization of this area, the four-gluon amplitude M^, to first 
order in y/XsYM, matches exactly the BDS ansatz, where the strong coupling scal- 
ing function is plugged in. This provides a remarkable check on the BDS ansatz as 
well as a hint towards the WL/scattering amplitude duality, since the strong coupling 
computation of the amplitude strikingly parallels that of a light-like WL. 

Motivated by the analogy with the four dimensional case and by evidence in favor of 
WL/amplitude duality and BDS exponentiation, we investigate the ABJM four-point 
amplitude at strong coupling, by following the same steps as in [38]. 

At strong coupling where the 't Hooft parameter A = N/K is large, and in the 
intermediate regime K N K 5 the AdS/CFT correspondence provides a dual 
description of the ABJM in terms of type IIA supergravity on AdS 4 x CP 3 . 

The dual background in string frame is 

ds 2 = — i^ds 2 AdS4 + ds 2 cp ,j (5.24) 

where in l s units the AdS± radius is given by R\dSi = {jc = ^ 25 ^ 2 R KN - = \/2 n \f\. 

A first indication that the general prescription for computing scattering amplitudes 
at strong coupling could still be (5.23) with the parameters conveniently adapted to 
the three dimensional model, comes from observing that the ratio of the two AdS radii 
coincides with the ratio of the scaling functions at leading order in the couplings. In 
fact, taking into account that for Af = 4 SYM at strong coupling 

f(XsYM) = + O (A°) (5.25) 
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whereas for the ABJM theory [19]-[30] 



fcs(\) = V2\ + O(\ ) 



(5.26) 



it is easy to see that 



R A 

R 2 A 



'■AdSk _ fi^SYAl) 
'■AdS 4 /cs(A) leading 



(5.27) 



Moreover, the string solution in AdS§ describing a four-point amplitude/light-like 
WL for Af = 4 SYM at strong coupling [38] may be straightforwardly embedded in 
AdS 4 as well. Therefore the four-point amplitude should be trivially readable from 
the Af = 4 SYM result (5.23) by changing the AdS radius. Indeed this supports the 
extension of the BDS-like ansatz (5.16) to strong coupling. 

The subtle point in identifying the Af = 4 solution with the ABJM one comes 
with regularization of infrared divergences. In [38], the strong coupling analogue of 
dimensional regularization is spelled out. This amounts to continuing the dimensions 
of the Dp-branes sourcing the AdS$ x S 5 background from p = 3 to p = (3 — 2e). 
Correspondingly, the new solution for the modified metric leads to the expression 



for the regularized world-sheet action that, once minimized, will provide the four- 
point amplitude (here c t is an e dependent constant and is the dimensionless 't 
Hooft coupling in dimensional regularization). 

In the context of ABJM we have not been able to find a similarly well-motivated 
regularization procedure 4 . However, guided by the analogy between the four-point 
ABJM and Af = 4 SYM amplitudes at weak coupling, we are tempted to employ the 
prescription (5.28) to regularize the action also in the AdS4 context. Following the 
same steps as for the Af = 4 SYM case, it leads to a strong coupling version of the 
ABJM four-point amplitude given by 




(5.28) 



M4 = e 



^+^(ln 2 (f) + ^)+Consts+o(^) 



(5.29) 



where the leading infrared divergence is 




(5.30) 



4 Although cutoff regularization works fine in three dimensions, we prefer to insist on a dimensional- 
like one in order to compare with our expression (4.54). 



Even though this prescription lacks strong motivations, it definitely captures the essen- 
tial features of the amplitude, such as the leading singularity and the coefficient of the 
non-trivial finite piece, which matches the strong coupling value of the ABJM scaling 
function (5.26). 

The generalization to the ABJ model is not straightforward. Here the situation 
is slightly subtler, since concerns have arisen on the integrability of the corresponding 
o"-model in the dual description [18]. Nevertheless, to first order at strong coupling 
we still expect the amplitude to be described by (5.29). The reason is that unitarity 
requires I — \M — N\ < K [2]. Hence at strong coupling, where M, N ^> K, the shift in 
the ranks is negligible compared to v^A of ABJM, and it affects the solution at higher 
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orders only, starting from O ( -73= ) [66] 
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A. Notations and conventions 

We work in three dimensional euclidean Af = 2 superspace described by coordinates 
(x»,6 a f), a, /3 = 1,2 . 

Spinorial indices are raised and lowered as (we follow conventions of [67]) 



^ = ^ = ^C Pa (A.l) 

where the C matrix 

**=(°;) c -=(?o) < a - 2 » 

obeys the relation 

C a ? C 1& = 8 a 1 5? s - S a s 5^ (A.3) 
Spinors are contracted according to 

^x = ^ a x a = x a *P a = xi> ij 2 = ^ a ij a (A.4) 

Dirac (^) a ^ matrices are defined to satisfy the algebra 

(tT 7 (tT/» = -gTS a p + i <r> fr p ) a p (A.5) 

Trace identities needed for loop calculations can be easily obtained from the above 
algebra 

tr(nn = (Yfa = -2<T (A.6) 

tr(rrr) =-m%(i v f,(i p )\ = ^^ p (a.7) 
= nsr<r-<r<r+<r<r') (a.s) 

Using these matrices, vectors and bispinors are exchanged according to 

coordinates: x M = (7 M ) Q( g x Q/3 = \ (7 M ) a/3 

derivatives : <9 M = \ (7^ d afi d a/3 = irf)afi d n ( A -9) 

fields: A^^^A^ Aa p = j= 2 (YUA, 

It follows that the scalar product of two vectors can be rewritten as 

p-k = \ V ^K^ (A.10) 
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Superspace covariant derivatives are defined as 

D a = d a + l -f d afj , D a = d a + l -6?d a p (A.n) 

and satisfy the anticommutator 

{D a ,D p } = id aP (A.12) 

The components of a chiral and an anti-chiral superfield, Z(xl,9) and Z(x R ,9), 
are a complex boson 0, a complex two-component fermion ip and a complex auxiliary 
scalar F. Their expansions are given by 

z = ( f ) (x L ) + e a ^ a (x L )-e 2 F(x L ) 

Z = 4>{x R ) + T$ a {x R ) - t F{x R ) (A.13) 

where x^ L = x^ + i9~f ti 9, x R — x^ — iQ^Q. 

The components of the real vector superfield V(x, 9, 9) in the Wess-Zumino gauge 
(V| = D a V\ = D 2 V\ = 0) are the gauge field A a p, a complex two-component fermion 
A Q , a real scalar a and an auxiliary scalar D, such that 

V = i 9 a 9 a a(x) + 9 a f V2A a p{x) - 9 2 T\ a (x) - f 9 a X a (x) + 9 2 1 D(x) (A. 14) 



The vector superfields (V, V") are in the adjoint representation of the two gauge groups 
U(M) x U(N), that is V = V A T A and V = V A f A , where T A are the U(M) generators 
and T A are the U(N) ones. 

The U(M) generators are defined as T A = (T°,T a ), where T° = ^= and T a 
(a = 1,...,M 2 — 1) are a set of M x M hermitian matrices. The generators are 
normalized as Tr(T A T B ) = 5 AB . The same conventions hold for the U(N) generators. 

For any value of the couplings, the action (2.1) is invariant under the following gauge 
transformations 

e v -)■ e iKl e v e- iM -)■ e^e^e^ (A.15) 



where A 1; A 2 are two chiral superfields parametrizing U(M) and U(N) gauge transfor- 
mations, respectively. Antichiral superfields transform according to the conjugate of 
(A. 16). The action is also invariant under the U(1) R R-symmetry group under which 
the A 1 and Bi fields have |-charge. 
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B. Integrals in dimensional regular izat ion 

In this Appendix we list a number of properties for momentum integrals entering the 
evaluation of four-point scattering amplitudes. 

We work in dimensional regularization, d = 3 — 2e, with dimensional reduction 
(spinors and e-tensors are kept strictly in three dimensions). 

B.l Tools 

When evaluating one-loop bubbles and two-loop factorized and composite bubbles we 
used the G[a, b] functions defined by 

1 T(a + b-d/2)T(d/2-a)T(d/2-b) 

GMJ "(4^ r(a)r(6)r(d-a-6) • 

Moreover, we introduce the compact notation r(ai|...|a„) = T(a 1 )...T(a n ). 

In order to write a momentum integral in its Feynman parametrized form, the basic 
identity is 

1 1 I>! + ••• + «„) } d(5 1 ...d(5 n 5([5 1 + --- + (5 n -l)ff 1 - 1 ...(5«-- 1 



A^'"A^ T( ai \ ... K) J ( /3l A 1 + ... + /3 n A0 ai+ - +a " 

o 

(B.2) 

where Aj are generic propagators. The integration over Feynman parameters makes 
often use of the identity 

i 

[dp,... d(3 n Sfa + ... + &- l) ^ . . . z^- 1 = iMm^L (B.3) 

J 1 («i H \-a n ) 

o 

The most complicated computations of the one and two-loop on-shell amplitude 
were performed using Mellin-Barnes representations [68, 69]. These representations are 
based on the identity 

inn 

,2 \ s 



f dsT(-s)T(s + a) (^-] , (B.4) 



(k 2 + M 2 ) a (M 2 ) 

where the contour is given by a straight line along the imaginary axis such that in- 
dentations are used if necessary in order to leave the series of poles s = 0, 1, • • • , n to 
the right of the contour and the series s = —a, — a — 1, • • • , — a — n to the left of the 
contour. 
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After Feynman-parametrizing a triangle integral and using (B.4), the following 
formula holds 

d d k 1 1 



(27r) d k 2 ^ (k - p) 2 ^ (k + q) 2 ^ (4vr) rf / 2 \[ T{m)T{d - J2i Hi 



■x 



r ds dt r (-s|-t|f -//i-// 2 -s|f -^i-/i 3 -t|^i+s+t|^/ii-f + s +t) 

X 7 (2ni) 2 WF* (q 2 )^ (P + g ) 2 ( B + t +EiW-«i/2) ^ • ^ 

— ioo 

while for a vector-like triangle we have 
d d k k u 



I 



(2n) d k 2 ^ (k - p) 2 ^ (k + q) 2 ^ 

ioo 

(4vr)- d / 2 f dsdt r(-s|-t|Afi+s+t|E^-|+s+t) 



x 



Ilir(A*i)r(d-X>i + i) 7 (2«) 2 ( p 2 )-s( ? 2)-t( p + ? )2(s + t + E iW -<i/2) 

—ioo 

[r(f-/i 1 -/i 2 -s|f-/i 1 -/i 3 -t + l)p 1 '-r(f-/i 1 - /U2 -s+l|f-/i 1 - /U 3-t)g"] (B.6) 

where the multiple contours are taken using the convention already mentioned for the 
relative position of the poles. When the position of a pole is chosen differently compared 
to the convention, it is customary to use the notation T*(z) for the gamma function 
involved. 

We can proceed along similar lines for writing the Mellin-Barnes representation 
for box diagrams, vector-like boxes, etc. Using these representations, along with the 
Barnes first lemma 



JOO 



ds ^ 1T ,„ 1TVJ , T(a + c)T(a + d)T(b + c)T(b + d) 

_ T (a + ,)r(6 + s)r(c - s)r(d - ,) = -> ' T \ a + h ' + \ + j 1 CB.7) 



/ 

—ioo 

we have been able to compute one and two-loop amplitudes in a manifestly analytical 
way, without performing numerical evaluations. 

B.2 Bubbles 

At one-loop, the evaluation of simple bubbles is required. Feynman parametrizing the 
integrand and working off-shell (p 2 ^ 0), we easily obtain 

/d 3 ~ 2e k 1 1 1 

W^wj^ = Gll ' 1] W*->~W\ +oie) (B ' 8) 

where we have defined \p\ = a/p 2 and G[a, b] is given in (B.l). On the other hand, if we 
are on-shell, p 2 = 0, the integral reduces to a tadpole-like integral and in dimensional 
regularization it vanishes. 
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B.3 Triangles 



We begin by evaluating the scalar triangle diagram of Fig. 5. 




When the external momenta are off-shell (p 2 7^ 0), the integral can be computed 
in three dimensions with no need for regularization. Since for D = 3 the triangle with 
propagator exponents (1, 1, 1) satisfies the uniqueness condition [70], it evaluates to a 
rational function 



d 3 k 



1 



(2tt) 3 p (k - Pl y (k + P2 ) 2 
1 



(B.9) 



The corresponding integral, when evaluated on-shell (p 2 = 0) and in dimensional reg- 
ularization, can be easily treated by Feynman parametrization and is given by 



T(pi,p 2 ) 



d 3 ~ 2e k 



(2tt) 3 - 2£ k 2 (k - pi) 2 (k + P2) 2 
r(3/2 + e)r 2 (-l/2 -e) 1 



(47r)f-T(-2e) 
1 + 0(e 2 ) 



bi+p 2 | 3+2e 



2|pi +P2I 



(B.10) 



Therefore, in dimensional regularization and on-shell limit, the scalar triangle can be 
set to zero. 



We then consider the on-shell, vector-like triangle. Again, by Feynman parame- 
trization, it is straightforward to show that the integral is given by 



7^ a/3 (pi,p 2 ) 



d^k 



k a/3 



(2vr) 3 - 2e k 2 (k - Pl ) 2 (k + p 2 ) 2 



p =o 



r(3/2 + e)r(l/2-e)r(-l/2-e) ( Pl -p 2 



,a/9 



e-5-0 



> B( Pl + p 2 ) 



(47r)3-T(l-2e) 
(P2 - 



(B.ll) 



where £>(pi + p 2 ) is the bubble in eq. (B.8). 

An important observation is that, as a consequence of the on-shell conditions, the 
vector-like triangle satisfies the following identities 



7v(Pl,P2) • (P1+P2) = -TvijPuPl) ■ {Pz+Pi) = 

B.4 Boxes 

We now consider scalar and vector-like box diagrams drawn in Fig. 6. 

Pi. .V± 



k-pi 




k +p 4 



k — Pl -p 2 

P2 P3 
Figure 6: The box diagram. 

In terms of the Mandelstam variables, the scalar integral is written as 

(P~ 2e k 1 



Q(M) 



(27r) 3 - 2e k 2 (k - Pl ) 2 (k — pi — p 2 ) 2 (k + p 4 ) 2 
while the vector-like integral is 

d 3 - 2e k k" 



Qv 



(2vr) 3 - 2e k 2 (k - pi) 2 {k - pi- p 2 ) 2 (k + Pi ) 2 



(B.12) 



(B.13) 



(B.14) 
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The scalar integral can be evaluated at leading order in e. Feynman parametrizing 
the integrand in eq. (B.13) we obtain 

r>( A - r ( 5 / 2 + e ) f fyi d Vz d Vz dy 4 5(Y Ji yi- 1 ) {n , ,\ 

J{S, ) ~ (4^)3/^7 ( yi y 3 s + y 2 y 4 t)^ ^ 

Expressing the denominator as a Mellin Barnes integral 

1 1 f du {yiVzs) 1 



( ^r(-«)r(« + 5/2 + e ) 7 {yim £+ (B.i6) 



(yiy 3 s + l/2|/4t) 5 / 2+e r(5/2 + e)y 2m 
and integrating on the Feynman parameters we obtain a one-fold representation 

2r(J +2f) 1 dU ^r(-u)r 2 (-3/2-e-u)r(5/2 + e+u)T 2 (l+ U )X u (B.17) 



(47r) D / 2 r(l -2e)t 5 / 2 + e J 2ni 

We note that the MB integral is multiplied by an e factor and the integral itself is well 
defined when e — > 0. Therefore, to leading order in e we have 

2M) = 4^572- / ^f{-u)T\^/2-u)T^/2 + u)T\l + u)X- + 0{e 2 ) 

^^^(fdX) + h(X)) + 0(e 2 ) (B.18) 

where X = s/t. By closing the contour on the right, fi(X) is the sum of the residues at 
the poles of F(—u), whereas fi{X) is the contribution from the poles of T 2 (— 3/2 — u). 
The fi(X) function is easily computed and gives 



aw - * 2 E t^> - ^ 3/2 (4P ^ " < B19 > 

The $2{X) function is more complicated since we have to deal with double poles. A 
first set of manipulations leads to 



n=0 



where ^(x) is the digamma function. 

The first term in (B.20) is easily summed to 



^(-X)T( l/2 + „) ^ 



n=0 n 
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The second series in (B.20) can be summed by using the trick 



(-X)"T(-l/2 + n)#(-l/2 + n) _ d_ I (-X) n T(a + n) 

^-^ n\ dn \ J—~< 



n=0 



K n=0 



o=-l/2 



d 



T(a) 



2^VTTX (ln(l + X) - #(-1/2)) 



da y(i + x)* /a= _ 1/2 

For the third term in (B.20) we use the following identity 



(B.22) 



#(l + n) = -j E 



1 -t n 
1-t 



dt 



^B.23) 



to rewrite the digamma function inside the series. By exchanging the order of the sum 
and the integral and summing the series, we obtain 



E 

n=0 



(-X) n r(-l/2 + n)#(l + n) 



2^ Vl + X - J dt 
-2vWTTX (#(-1/2) + 2 ln(l + 7T ^)) + 4v^F 



VTTx -VT+tx 



l-t 



(B.24) 



where the integral in the first line has been performed using Mathematica. 

Summing (B.21, B.22, B.24), after many non-trivial cancelations and simplifica- 
tions we obtain as a final result 



Q(s,t) 



(Vst) 



+ C(e 2 ) (B.25) 



The vector-like box integral (B.14) can be computed by using the same methods 
as before or, equivalently, by using Passarino-Veltman reduction to write it as a linear 
combination of scalar integrals. In any case, at leading order in e, we obtain 



Qv 



1 



In 



s + Vi+ Vs + t 



1 2( v / st) 3 \ VsTt x " \y/s + y/i- v^Tt 



s (pi - paT +t( Pl + p 2 y] 



(B.26) 



It is interesting to note that the projections of <2y i n the directions of p\ and p 4 become 
very simple 

e / 1 



Q -Pi = I (^2 " W2) + °^ 2 ) and Q. PA = l^-^ 2 



(B.27) 
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since the logarithm term drops. 

In fact, these two projections, can be calculated to all orders in e. Writing k-p% and 
k ■ j>4 in the numerator of (B.14) as the difference of two squares, the integral reduces 
to the difference of two triangles. Therefore, using the results of Subsection B.2, we 
obtain 

= r(3/2 + 6)r 2 (-l/2-e) ( 1 1 

-<_- 'if 1 



2 (47r)^/ 2 r(-2e) V* 3/2+e * 3/2+e 
_ r(3/2 + e)r 2 (-l/2- e ) / 1 1 \ 

- 2(4vr)^r(-2e) ^ / ' 1 J 

Since the leading term for e — >• coincides with the expressions (B.27), this is a consis- 
tency check of our results. 



C. The basis of dual conformally invariant integrals 

In this Appendix we give the proof that the linear combination Ic )S = I\ s — 1^ of scalar 
integrals defined in (5.1, 5.4), when evaluated in three dimensions, is free from IR 
divergences. This allows to conclude that the actual basis for two-loop amplitudes is 
lis, 12s, hs, hs plus their t-counterparts. 

We consider I is in (5.4) and apply the following identity 

X m = X 56 + X 45 + 2 ^45 • x 56 (CI) 

to its numerator, thus decomposing it as (we neglect the factor xf 3 ) 




where arrows indicate contractions of the corresponding variables at the numerator (see 
Fig. 3 for the labeling of momenta and dual variables). Here we already recognize the 
emergence of the infrared divergent integral I\ s , a triangle-box which could diverge, 
having unprotected cubic vertices, and a double-box whose cubic vertices are mitigated 
by the presence of a non-trivial numerator. 

We focus on the triangle-box and handle it by using the identity 

x 25 = x 56 + x 26 ~ 2^56 ' ^26 (C-2) 

in the numerator factor. The final result can be cast into the following graphical relation 
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where the IR divergence has been completely isolated in the last term. Therefore, taking 
the linear combination J ls — I 4s the divergences cancel at the level of the integrands 
and we are left with a well-defined dual conformally invariant integral. 
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